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The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the invariant dif- 
ferential cross section for production of Kg, uj, rj , and <j> mesons in p+p collisions at y/s = 200 GeV. 
Measurements of uj and 4> production in different decay channels give consistent results. New results 
for the uj are in agreement with previously published data and extend the measured pr coverage. 
The spectral shapes of all hadron transverse momentum distributions measured by PHENIX are well 
described by a Tsallis distribution functional form with only two parameters, n and T, determining 
the high-pT and characterizing the low-pr regions of the spectra, respectively. The values of these 
parameters are very similar for all analyzed meson spectra, but with a lower parameter T extracted 
for protons. The integrated invariant cross sections calculated from the fitted distributions are found 
to be consistent with existing measurements and with statistical model predictions. 



PACS numbers: 25.75.Dw 

I. INTRODUCTION 

The PHENIX experiment at the Relativistic Heavy Ion 
Collider (RHIC) at Brookhaven National Laboratory has 
measured the production of a wide variety of hadrons (7r, 
K, r], rf , uj, <f>, p, J/ip, and ip') at midrapidity in p + p 
collisions at y/s = 200 GeV. The measurements were 
performed using a time-of-flight technique for charged 
hadron identification and via reconstruction of various 
photonic, hadronic, and dielectron decay modes for neu- 
tral hadrons. The measured transverse momentum spec- 
tra extend over the range from zero to 20 GeV/c. Precise 
measurements of hadron production in p + p collisions 
are crucial for a deeper understanding of QCD phenom- 
ena such as parton dynamics and hadronization. They 
also provide a valuable baseline for particle and jet pro- 
duction in heavy ion collisions, essential to the needs of 
the RHIC heavy ion program. 

There exists a large body of experimental data on 
hadron production in p + p collisions measured at the 
ISR, SppS, Tevatron, and RHIC [Ml. At high p T the 
spectra display a power law behavior that becomes more 
and more evident as the interaction energy increases. In 
this regime, the spectra are well described by pertur- 
bative QCD together with measured proton structure 
functions [20j. At low pr, typically pr < 2 GeV/c, a 
region which accounts for the bulk of the produced par- 
ticles, the spectra are governed by processes that belong 
to the non-perturbative regime of QCD and are not yet 
fully understood. In this px region, the spectra reveal 
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an exponential behavior which can be explained with the 
assumption that secondary particles are emitted from a 
thermalized system with at most short-range correlations 
and obeying Boltzmann-Gibbs statistics [2l|. In this ap- 
proach, the inverse slope parameter T can be interpreted 
as the temperature of the system. However, that would 
require some mechanism of local thermal equilibrium in 
p + p collisions which is not yet established. It is also 
known that the particle spectra are best described by an 
exponential in mj- rather than in px [22J . According to 
the observation that the temperature parameter T is the 
same for different particles a spectral shape is also the 
same, the so-called mr-scaling observation |18l l23j. 

The two regimes described here, and the pr-region 
where their contributions are predominant are commonly 
designated as "soft" and "hard" . There is no obvious 
boundary between them and the question as to what ex- 
tent the mechanisms in each region are distinct is difficult 
to address. In this paper it is shown that the spectral 
shapes of all hadrons produced in p + p collisions at ^fs 
=200 GeV measured by PHENIX are well described by 
one single distribution without making a distinct divi- 
sion into two regions. The Tsallis j24j distribution, also 
referred to as Levy distribution fA |25|. has only two pa- 
rameters, T and n, that characterize the low- and high-p^ 
regions of the spectra, respectively. This distribution has 
been shown by Tsallis to result from a postulated gen- 
eralization of the Boltzmann-Gibbs entropy. It has been 
suggested to be relevant for various types of systems, such 
as those with long-range correlations, or non-ergodic fill- 
ing of the available phase space. Boltzmann statistics 
and exponential distributions are recovered in the limit 
that correlations disappear. The parameter T then re- 
covers the usual interpretation as the temperature of the 
system. 



4 



In a number recent publications the Tsallis statisti- 
cal distribution was successfully applied to describe the 
Heavy Ion and High Energy data in over a wide range of 
incident energies and centralities (26l - [30| . Physical mech- 
anisms responsible for the successful application of the 
non-extensive statistical approach to the description of 
the particle spectra in this systems is a topic of discus- 
sion [31H36J. The analysis presented in this paper uses 
the Tsallis formalism as a parameterization to describe 
the particle spectra and compares it with other parame- 
tcrizations used for the spectra approximation. Common 
features and differences revealed in such an approach are 
data driven and should contribute a better understanding 
of particle production mechanisms. 

The successful description of the particle spectra with 
the Tsallis distribution allows to accurately calculate the 
integrated particle yield and mean momentum, even for 
species measured only in a limited momentum range. 
The integrated particle abundances provide important 
information on the bulk properties of the soft particle 
production. In particular, the comparison of the parti- 
cle yields to statistical model predictions can be used to 
infer the degree of chemical equilibration. In the case 
of heavy ion collisions, the success of statistical model 
fits to the particle yields [13, HI] suggests that chemical 
equilibration is essentially complete. These models have 
also been used to describe particle production in p + p 
collisions |!! S(J • 

In this paper we present new PHENIX results on the 
production of neutral mesons in p + p collision at y/s 
=200 GeV and compare the PHENIX data with the pa- 
rameterizations commonly used to describe particle spec- 
tra in relativistic p + p collisions, including the Tsallis 
parametrization. It is demonstrated that the latter one 
describes the data in the entire range of measured pt 
most accurately. The parameter values extracted from 
the fits are given for all measured particles. 

The paper is organized as follows: Section |TT] gives 
a description of the PHENIX experimental setup and 
detector subsystems. Section Mil describes the analy- 
sis methods used to measure the transverse momentum 
spectra of different hadrons for p + p collisions at y/s 
=200 GeV. In Section HVl the properties of the measured 
transverse momentum spectra are analyzed. In SectionlVl 
the scaling properties of the particle spectra are discussed 
and the calculated integrated yields are compared with 
published results and with statistical model calculations. 
The measured invariant cross sections are tabulated in 
tables given in the Appendix. 



II. PHENIX DETECTOR 

The PHENIX detector is designed as a high rate and 
fine granularity apparatus that utilizes a variety of de- 
tector technologies to measure global characteristics of 
the events, and to measure leptons, hadrons, and pho- 
tons over a wide range of transverse momenta. The ex- 



perimental setup consists of two central arm spectrome- 
ters each covering A0 = tt/2 in azimuth at midrapid- 
ity \rj\ < 0.35; two forward muon spectrometers with 
full azimuthal coverage in the rapidity interval 1.2 < 
\r}\ < 2.4(2.2) for the North (South) arm and a sys- 
tem of "global" detectors. Each spectrometer provides 
very good momentum and spatial resolution and parti- 
cle identification capabilities. The detailed description 
of the detector can be found elsewhere [4l[. The ex- 
perimental results presented in this paper were obtained 
using the central spectrometers and global detectors of 
the PHENIX experiment schematically shown in Fig. [TJ 

Reconstruction of charged particle tracks and momen- 
tum measurements are performed with the drift cham- 
bers (DC) and the first layer of the pad chambers (PCI). 
The fiducial volume of the DC is located outside of the 
analyzing magnetic field of the detector and has an inner 
radius of 2.02 m with an outer radius of 2.46 m. Multiple 
layers of wires measure the track position with an angular 
resolution of ~ 0.8 mrad in the bending plane perpendic- 
ular to the beam axis. The PCI, located just outside the 
outer radius of the DC, has a spatial resolution of a$ ~ 
2.4 mm and a z ~ 1.7 mm and provides the z-coordinate 
of the track at the exit of the DC. The momentum of 
a particle is determined by the measured bending an- 
gle in the axial magnetic field of the central magnet [42[ 
assuming that the particle originates from the collision 
vertex. The DC momentum resolution is estimated to be 
Sp/p = 0.7 © l.l%p [GeV/c]. Track matching with hits 
in the second (PC2) and third (PC3) pad chamber layers 
located at radii of 4.2 m and 5.0 m, respectively, rejects 
tracks from secondaries originating either from decays of 
long-lived hadrons, or from interactions with the struc- 
ture of the detector. Such tracks have not passed through 
the full magnetic field and therefore have improperly de- 
termined momenta that is typically overestimated. A 
detailed description of the PHENIX tracking system can 
be found in 0,S3- 

The primary purpose of the PHENIX Electromagnetic 
Calorimeter (EMCal) is to measure the position and en- 
ergy of photons and electrons. The EMCal covers the full 
acceptance of the central spectrometers and is divided 
into eight sectors. Six of the EMCal sectors located at 
the radius of 5.0 m are built of lead-scintillator (PbSc) 
and consist of 15552 individual towers with a granularity 
of 5.5x5.5 cm 2 and a depth of 18 Xq. The two other sec- 
tors located at the radius of 5.2 m are built of lead-glass 
(PbGl) and consist of 9216 lead-glass Cerenkov towers 
with a granularity of 4x4 cm 2 and a depth of 14.4 Xq. 
Due the fine segmentation of the EMCal the electro- 
magnetic showers typically spread over several towers. 
This spread provides the means to analyze the posi- 
tion and shape of the shower, and to reject hadrons 
which produce showers of a different shape. The spatial 
resolution of the PbSc(PbGl) EMCal sector is a(E) = 
1.55(0.2) © 5.7(8.4)/ y/E[GcV] mm for particles at nor- 
mal incidence. The energy resolution of the PbSc(PbGl) 
calorimeter is SE/E = 2.1(0.8)% © 8.1(5.9)/ y/E[GeV}%. 
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West Beam View East 



FIG. 1: Schematic view of the PHENIX central spectrometers and particle decay modes analyzed in this paper. 



The Time-of-Flight (TOF) subsystem is used for 
hadron identification based on momentum measurements 
in the DC and PCI combined with flight path length 
from the collision vertex (45[. The TOF is located be- 
tween the PC3 and the PbGl at the radius of 5.0 m and 
covers about 1/3 of the acceptance of one central arm. 
The TOF detector consists of 10 panels each contain- 
ing 96 segments equipped with plastic scintillators and 
photomultiplier readout from both ends. The time res- 
olution of ~ 120 ps enables ir/K and K/p separation in 
the transverse momentum range 0.3-2.5 GeV/c and 0.3- 
5.0 GeV/c, respectively. The lower limit is denned by the 
energy loss of different particles in the detector material. 

The Ring- Imaging Cerenkov (RICH) is the primary de- 
tector for e/7r separation. It provides an e/ir rejection 
factor of ~ 10 -3 for tracks with momenta below the pion 
Cerenkov threshold of ~ 4 GeV/c in the CO2 used as a 
radiator gas. The RICH detector in each arm has a mir- 
ror measuring 20 m 2 that focuses the light onto an array 
of 2560 photomultipliers. The material of the PHENIX 
central arm that precedes the RICH has been kept to 
just ~ 2% of a radiation length in order to minimize the 
background contribution of electrons from 7-conversion. 
The PHENIX RICH and TOF detectors are described in 
more detail in [451 ] . 

The Beam-Beam Counters (BBC) are used for trig- 
gering, determination of the collision time, and location 
of the vertex along the beam axis, z v tx [46j- The BBC 



consists of two sets of 64 Cerenkov counters surrounding 
the beam pipe, and located at a distance of ± 1.44 m 
from the center of the interaction region. Each counter 
covers the full azimuth and the pseudorapidity interval 
3.1 < \r}\ < 3.9. The z-coordinate of the collision vertex 
is determined with a typical resolution of 2 cm in p + p 
collisions by the timing difference of the signals measured 
by the two sets of BBC counters. The time average of 
all BBC counters gives a start time for the time-of-flight 
measurements. The minimum bias trigger in p + p colli- 
sions is generated when at least one counter fires in each 
BBC set of counters, and the collision vertex calculated 
on-line is |z v tx| < 38 cm. The efficiency of the minimum 
bias trigger is estimated to be (55±5)% of the total in- 
elastic cross section of <J^ el = 42 ± 3 mb. Further details 
about the BBC subsystem of the PHENIX detector can 
be found in [46j. 

Due to the high rate of p + p collisions at RHIC 
PHENIX employs several specialized triggers which en- 
able the experiment to sample more of the delivered lumi- 
nosity for rare events. Besides the minimum bias trigger, 
the experimental results presented in this paper were ob- 
tained using the EMCal-RICH Trigger (ERT). 

The EMCal is used to trigger on rare events with large 
energy deposit originating primarily from high energy 
photons or electrons. The analog sum of signals from 
4x4 adjacent towers is compared with a trigger thresh- 
old of 1.4 GeV. In addition, a combination of the EMCal 
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and the RICH signals is used to build the ERT trigger 
which is designed to select events containing electrons. 
The trigger fires when the analog sum of signals from 
2x2 adjacent towers in the EMCal exceeds a threshold of 
0.4 GeV (setting used in the 2005 physics run) or 0.6 GeV 
(used in 2006) in geometrical coincidence with a signal in 
the associated RICH trigger tile (4x5 PMTs) determined 
using a look-up table. 



III. NEUTRAL MESON MEASUREMENTS 



In this section we describe the analysis details of the 



Kg — > 7r"7r", lu — > 7r"7r T 7r~, lu — > 7t"7, lu — > e T e~, 
r/' — > r)ir + TT~, <f> — > K + K~, and <j) — > e + e~ measure- 
ments in p + p collisions at ^/s —200 GeV. These mea- 
surements complete and extend previous neutral meson 
spectra results measured by the PHENIX experiment and 

published in H S i & Ea Ei 

The measurements are based on a data sample repre- 
senting a total integrated luminosity of 2.5 pb^ 1 within a 
vertex cut of |z v tx| < 30 cm accumulated by the PHENIX 
experiment in 2005. The data were collected using mini- 
mum bias and ERT triggers. 




Invariant Mass (GeV/c ) 



A. Reconstruction of neutral mesons 

Here we discuss the analysis details and main parame- 
ters of the invariant mass distributions reconstructed for 
different decay modes. 



1. Selection of the n , r\ — > 77 candidates 

Most particles studied in this section decay producing 
a 7T° or r\ meson in the final state, which in turn decays 
into a 77 pair at the point of primary decay. The analysis 
procedures for the measurement of the inclusive ir° and 
77 invariant transverse momentum spectra in p + p colli- 
sions have been published previously 0, H, Elf . Me- 
son candidates were reconstructed from pairs of clusters 
in the EMCal with energy E~ t > 0.2 GeV assuming that 
they correspond to photons originating from the collision 
vertex. A shower profile cut was used to reject broader 
showers predominantly produced by hadrons 49]. The 
invariant mass distribution for cluster pairs is shown in 

Fig.m 

The width of the peaks is determined largely by 
the EMCal energy resolution. For 7r°(?7) meson can- 
didates the width decreases from 12(40) MeV/c 2 to 
9(30) MeV/c 2 between 1 GeV/c and 3 GeV/c of the pair 
transverse momentum. 

The reconstructed positions and widths of the peaks 
are in agreement with simulation results once detector 
resolution and trigger biases have been taken into ac- 
count. The measured mass peaks were parameterized as 



FIG. 2: Invariant mass distribution for 77 pairs in the range 
4 < pr (GeV/c) < 6. The inset shows an enlargement of the 
region around the r\ mass. 



a function of the 77 pair pt- For further analyses involv- 
ing 7T° or r\ mesons in the final state we selected pairs 
with pt > 1 GeV/c and an invariant mass within two 
standard deviations of the measured peak position. All 
77 pairs satisfying these criteria were assigned the nomi- 
nal mass of the meson and the photon energies were 
rescaled by the ratio of the nominal to the reconstructed 
masses. 



2. 7r°7 and Ag -> ttV 

The reconstruction of lu — » 7r°7 and K$ — > tt°tt° decays 
was performed by combining ir° candidates with either 
all other photons with energy > 1 GeV \4\ or with 
all other ir° candidates from the same event. Combina- 
tions using the same EMCal clusters more than once were 
rejected. 

Invariant mass distributions for 7r°7 and ir°ir° decays 
are shown in Fig. [3] The width of the 10 meson peak is 
~ 30 MeV and has a weak pt dependence. The width 
of the Kg peak is ~ 15 MeV. The signal-to-background 
ratio (S:B) increases from 1:30 (1:4) to 1:5 (1:2) for ui 
(K s ) mesons as the transverse momentum increases from 
2 to 12 GeV/c. 

The main difference in the analysis of the co and Kg 
decays was due to the large lifetime of the A'^-meson. 
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FIG. 3: Invariant mass distribution for 7r°7 (left) 
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and 7T 7r° (right) decays at 4 < p T (GeV/c) < 6. 



Neutral pions coming from the decays of high-pT K s 
originate from a displaced vertex and their reconstructed 
mass and width need to be parameterized in a different 
way compared to pions coming from the primary event 
vertex. This effect was studied using the PHENIX Monte 
Carlo. The correction was based on the mass and width 
of 7T° 's coming from kaon decays with a realistic pt dis- 
tribution, and on 7r°'s produced at the collision vertex 
with the inclusive pr distribution. 

3. LO,T] — ?> 7T°7r + 7r~, 77' — > rjTY + TY~ 

For the reconstruction of to, rj — > 7r°7r + 7r~ and 7/ — » 
7/7r + 7r~ decay modes we combined 7r°(?/) candidates 
with all pairs of oppositely charged tracks in the same 
event [1, 0]. Charged tracks accepted for this analy- 
sis were required to have momenta in the range 0.2 < 
Pt (GeV/c) < 8, and were assigned the charged pion 
mass. Tracks with momentum below 0.2 GeV/c do not 
go through the entire detector due to their large bending 
angle in the axial magnetic field of the central magnet. 
Tracks that appear to have momenta above 8 GeV/c are, 
for the most part, low momentum secondaries coming 
from the decay of long lived primaries. Because they do 
not originate from the collision vertex, their momenta are 
not calculated correctly. Invariant mass distributions for 
ir° (r])Tr + ir~ triples are shown in Fig. [4] The two peaks 
in the distribution shown in the left panel of the figure 
correspond to decays of 77 and w mesons. The width of 
~ 8 MeV/c 2 for the reconstructed r\ meson peak is similar 



to that of the rf meson peak shown in the right panel of 
Fig. The width of the w-meson peak is ~ 17 MeV/c 2 
which is narrower than that in the ui — » ir°j decay mode. 
This is due to the smaller difference between the masses 
of the primary particle and their decay products and to 
the better momentum resolution of the tracking system 
as compared to the EMCal in this momentum range. The 
signal-to-background ratio in the range of measurements 
changes from 1:10 (1:5) to 1:3 (1:2) for u> (r) 1 ) mesons. 
More details on the analysis of r\ and w-mesons can be 
found in [Hi]. 

Reconstruction of the <f> — > K + K~ decay was done 
by combining pairs of oppositely charged tracks. The 
tracks were required to have a momentum in the range 
0.3 < pt (GeV/c) < 8. Each track was assigned the 
charged kaon mass. Invariant mass distributions were 
accumulated in two different configurations: i) combining 
all tracks reconstructed in the PHENIX tracking system; 
ii) combining all tracks of one sign with tracks of the 
opposite sign identified as a kaon in the TOF subsystem. 
Examples of the invariant mass distributions for the two 
cases are shown in the left and right panels of Fig. [5l 
respectively. 

The use of particle identification improved the signal- 
to-background ratio by a factor of more than two at the 
expense of a more limited acceptance, resulting in a fac- 
tor of five loss in statistics. At low and intermediate pt. 
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FIG. 4: Invariant mass distribution for 7r°7r + 7r (left) and t]-k + -k (right) triplets in the momentum range 4 < px (GeV/c) < 6. 
The insert shows the invariant mass distribution after the background removal explained in sec. IIII Bl 



where the combinatorial background is high but the data 
sample has large statistics this method is preferable. The 
method without particle identification was more effective 
at intermediate and high px because of the significant 
gain in the acceptance. The highest px reachable with 
this method is limited by the available statistics in the 
minimum bias data sample. The two methods described 
here use different detector subsystems and produce dif- 
ferent shapes of combinatorial background and signal- 
to-background ratios. Use of the two methods allowed 
to extend the px coverage of the measurement and pro- 
vided a consistency check between the results obtained 
in the overlap region between 1.5 GeV/c and 4.5 GeV/c. 
The signal-to-background ratio changes from 1:10 to 2:1 
depending on the analysis method and the px bin. More 
details on this measurement can be found in [51| . 



5. cj, (j> — ^ e e 

Electrons are reliably identified by the PHENIX de- 
tector in the momentum range 0.2 < px (GeV/c) < 4. 
Electron identification is accomplished using the informa- 
tion from the RICH and EMCal subsystems by requir- 
ing at least two RICH phototubes to fire within the ring 
shaped area associated with a charged track. In addi- 
tion, the ratio of the associated cluster energy measured 
in the EMCal to the momentum measured in the track- 
ing system must satisfy \E/p — 1| < 0.5. The invariant 
mass distribution obtained by combining identified e + 



and e~ pairs is shown in Fig. [5] for pairs in the range 
0.5 < px (GeV/c) < 0.75. The two peaks correspond 
to ui + p and <fi mesons. The widths of the u> (</>) me- 
son peaks vary from 6.1(6.0) MeV/c 2 to 9.0(11) MeV/c 2 
from the lowest px to the highest px of the electron pairs. 
The signal-to-background ratio in the region of the uj (jfi) 
meson peaks changes from 1:2(2:1) to 3:1(6:1). 



B. Raw yield extraction 

To extract the raw yields the invariant mass distribu- 
tions near each peak were parameterized as the sum of 
signal and background contributions. 

For the signal, we used a Breit-Wigner function con- 
volved with a Gaussian function (BW*G). The Breit- 
Wigner describes the natural shape of the measured res- 
onance and the Gaussian takes into account the detec- 
tor resolution. Depending on the decay channel being 
analyzed, one or the other contribution may dominate, 
e.g. the Gaussian part is more important in decays like 
uj — > 7r°7, or — > 7r°7r°, and the Breit-Wigner part in 
decays like <j> K + K~ or ui, <fi —> e + e~ . 

In most cases the parameters of the BW*G function 
when fitted to the data were consistent with the values 
expected from simulation. In the highest px bins, where 
the available statistics becomes a limiting factor, we con- 
strained the Gaussian width based on simulations. 

The (j> —> K + K~ decay mode was treated somewhat 
differently. Kaons decaying in flight before passing com- 
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FIG. 5: Invariant mass distribution for K + K~ accumulated without particle identification (left) and with one 7f-meson 
identification (right) in the momentum range 4 < pr (GeV/c) < 6. The insert shows the invariant mass distribution after the 
background removal explained in sec. IIII Bl 



pletely through the PHENIX tracking system modify the 
shape of the invariant mass distribution compared to 
those passing through the detector without decays. This 
results in non-Gaussian tails of the detector response 
function, and thus the Breit-Wigner and Gaussian width 
parameters in the BW*G convolution mix together. To 
account for this effect a Monte Carlo sample was pro- 
duced with the natural width of the 4> set to zero and the 
kaon lifetime set to infinity. Using these samples allowed 
to disentangle the effects related to the kaon decays in 
flight. 

In the analysis it was verified that the peak positions 
and widths obtained from the fits to the data were in 
agreement with the simulated data to within 10%. In 
the measurement of the w, <j) — ¥ e + e~ decays other terms 
were added to the BW*G shape to account for p decays 
and for internal conversions taken from [53 . [53| . The 
contribution of p underneath the u peak was estimated 
using Breit-Wigner parameterization, with the assump- 
tion that the production ratio of p and w is 1 and in the 
fit their ratio was determined by their e + e~ branching 
ratios in vacuum equal to 1.53. 

To properly estimate the background under the peak it 
is necessary to assume that the shape of the background 
does not change rapidly. With this assumption one can 
expand the background shape in a Taylor series around 
the peak position and take the most significant terms of 
the expansion. A natural choice is to use a second or- 
der polynomial. The regions outside the resonance peak, 



where the background dominates, define the parameters 
of the fit. For a second order polynomial fit the back- 
ground varies smoothly under the peak. This may not 
be the case for higher order polynomial fits to the back- 
ground. 

The combinatorial background in the data has two 
main contributions. The first comes from the random as- 
sociation of uncorrelated tracks. Its shape is defined by 
the detector acceptance and the pt distribution of par- 
ticles in the event. This part of the background remains 
smooth in the mass interval comparable to the width 
of the peaks shown in Figs. [31151 The correlated part of 
the combinatorial background comes from partially or in- 
correctly reconstructed decays of true particles and jets, 
and may have a faster changing shape. In several anal- 
yses the most significant contributions to the correlated 
background were studied to verify that they do not affect 
the raw yield extraction procedure. For example, the de- 
cay rj — > 77 produces an 777 peak at around 0.6 GeV/c 2 
in the invariant mass distribution of 7r°7. Also, the de- 
cay of K® — > 7r + 7r~ produces a peak at ^1.07 GeV/c 2 
in the K + K~ invariant mass distribution when two pi- 
ons are erroneously assigned the kaon mass. In some 
cases these processes limit the mass range available for 
the background determination. The mass range used for 
determination of the background did not include regions 
where one could expect appearance of such peaks. 

The raw yields were measured in the following way. 
First, the invariant mass distributions in different px 
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FIG. 6: Invariant mass distribution for e + e— pairs in the 
momentum range 0.5 < pr (GeV/c) < 0.75. The insert shows 
the invariant mass distribution after the background removal 
explained in sec. MI Bl 



bins were fitted with the B W* G plus background in the 
mass range of ±5 combined widths of the BW*G around 
the nominal mass of the meson. The exact range var- 
ied slightly depending on particle species and the px 
bin. The background contribution, estimated by the sec- 
ond order part of the fit function, was subtracted from 
the measured invariant mass distribution and the result- 
ing histogram was used to count the raw yield. Bins 
lying within ±2.5 combined widths of the BW*G func- 
tion around the mass peak contributed to the yield. The 
same procedure was used to calculate the raw yield in 
the Monte-Carlo used for the acceptance evaluation. 

The systematic uncertainty of the raw yield extraction 
was usually the main contributor to the total systematic 
uncertainty. We evaluated this uncertainty by modify- 
ing the analysis procedure. The main goal of that was 
to change the shape of the background around the reso- 
nance peak in a manner similar to that shown in Fig. [5] 
To achieve this goal, analyses of the same decay modes 
were performed in different ways. For example, by re- 
quiring PC3 or EMCal hit matching for charged tracks, 
varying the minimum energy of 7 clusters, or modifying 
the selection criteria for 7r°(?7) candidates. Independent 
of that we also varied the parameters of the fit functions, 
such as the fit range and the order of the polynomial. 
Typically, six to ten raw yield values were accumulated 



for each px bin. After fully correcting each of them for 
the corresponding reconstruction efficiency the variance 
of the results was taken as the systematic uncertainty. 



C. Invariant mass resolution 

The invariant mass resolution of the detector plays an 
important role in the analyses described in this section. 
It depends on several factors. Use of the detector track- 
ing system or EMCal makes a large difference. The mo- 
mentum range of the analyzed particles is less important. 
The difference between the mass of the particle and its 
decay products contributes directly to the invariant mass 
resolution. To demonstrate this we consider the limiting 
case of a particle decaying into two massless products. In 
this case, one can approximate the invariant mass reso- 
lution with the simple relation 5m/ m = {l/\/2)8px/pT- 
The single particle momentum resolution was discussed 
in section [XT] Figure [7] compares this approximation with 
the widths of the peaks shown in Figs. -[5] The mea- 
sured widths are plotted as a function of the mass differ- 
ence between the particle and its decay products. The 
two lines in the plot are calculated for two body decays 
reconstructed either with the tracking system only, or 
with the EMCal only at pair px of 4 GeV/c. 
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FIG. 7: Invariant mass resolution of the PHENIX detector for 
different decay modes measured in the momentum range 4 < 
Pt (GeV/c) < 6, except for the e + e~ mode that is measured 
in the range 0.5 < pr (GeV/c) < 0.75. The lines indicate the 
expected detector mass resolution. 

As can be seen, the simple approximation describes the 
measured mass widths for the two body decays reason- 
ably well. The widths of the e + e~ decay modes are some- 
what narrower due to use of a lower momentum range. 
The results for the J/ip and ip' which are not shown in 
the plot are also consistent with the trend of the "track- 
ing" line. The <f> —> K + K~ represents the case where the 
assumption of massless products is least valid, neverthe- 
less, the agreement is still reasonable. 

The widths of the invariant mass peaks reconstructed 
with both the EMCal and the tracking systems are dom- 
inated by the EMCal resolution. However, due to the 
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energy correction applied to the 7 clusters forming ir° or 
77 candidates the widths of the peaks reconstructed with 
3 and 4 particles are below the "EMCal" line. 



D. Detector acceptance and efficiency 

1. Geometrical acceptance and the analysis cuts 

The determination of the detector acceptance was done 
using a single particle Monte Carlo simulation. Particles 
were uniformly generated within \y\ < 0.5 in rapidity 
and in full azimuthal angle. The range of the transverse 
momentum distributions were chosen to produce suffi- 
cient statistics in all pt bins for which the signal could 
be extracted from the data. For the acceptance calcula- 
tion the generated spectra were weighted to match the 
measured particle spectra. This procedure was done it- 
eratively. Kinematics of the three-body decays of the 77, 
u>, and 7/ mesons assumed the experimentally measured 
phase space density distributions |54l - [59| . 

A GEANT based simulation of the PHENIX detec- 
tor was tuned to reproduce the response of all detector 
subsystems and inactive areas. It was verified that the 
simulated positions and widths of the 7r°, K®, 77, u, 77', 
and <f> peaks were consistent with the values measured in 
real data at all pr's. The same analysis code was used 
for the reconstruction and analysis of the simulated and 
real data. 

The detector acceptance, calculated as the ratio of the 
number of fully reconstructed particles to the number of 
generated particles, is shown in Fig. [5J All curves take 
into account the detector geometry, particle decay kine- 
matics, performance of the detector subsystems includ- 
ing particle identification, and the analysis cuts. The 
efficiencies strongly depend on the particle momentum 
and rapidly decrease at low px for all species studied in 
this analysis, establishing a low px edge for the measure- 
ments. 



2. ERT trigger efficiency 

The analysis of several decay modes was based on data 
samples accumulated with the ERT trigger described in 
section [TTJ The ERT trigger efficiency was extracted us- 
ing the minimum bias event sample. Each EMCal cluster 
which set the ERT trigger bit to indicate a 7-cluster or 
electron was identified. The track or cluster had to also 
satisfy the analysis cuts of a particular decay mode, and 
match the region where the trigger bit was generated. 
The trigger efficiency was calculated as the energy spec- 
tra of such clusters divided by the energy spectra of all 
accepted clusters or electrons. Trigger efficiencies of pho- 
tons and electrons measured for one of the PbSc sectors 
as a function of cluster energy are shown in the left panel 
of Fig. E 



The trigger efficiencies grow steeply with energy, reach- 
ing 50% at values approximately corresponding to the 
on-line trigger threshold setting of 0.6 GeV for electrons 
and 1.4 GeV for photons. The curves saturate at approx- 
imately twice the threshold energy. The level of satura- 
tion is below 100% because of inactive areas of the ERT 
and the RICH efficiency. 

For the analyzed decay modes the trigger efficiency 
evaluation was done using the same Monte-Carlo sample 
as was used for the acceptance calculation. First we re- 
quired the particle to be reconstructed in PHENIX with- 
out ERT trigger requirement. Then, for all EMCal clus- 
ters associated to photons or electrons in the final state of 
the decay, we generated a random number between and 
1 and compared it to the magnitude of the curve shown 
in the left panel of Fig.[9]at the energy of the cluster. The 
particle was considered to fire the ERT trigger if at least 
one of the randomly generated numbers was lower than 
the corresponding value of the curve. The probability to 
fire the ERT trigger for all analyzed mesons is shown in 
the right panel of the same figure. 



3. Electron identification efficiency 

The electron identification efficiency is included in the 
acceptance efficiencies shown in Fig. [8] It was evaluated 
using a full detector Monte Carlo simulation which was 
tuned to adequately reproduce the RICH and the EMCal 
detector responses. To ensure that the electron identifi- 
cation efficiency was properly done in the simulation it 
was confirmed to agree with the efficiency measured with 
real data. 

For this comparison the data samples accumulated 
during special PHENIX runs were used. In those runs 
a 1.7% radiation lengths brass converter was installed 
around the RHIC beam pipe in the PHENIX interac- 
tion region. In this sample we selected electrons of 
both signs using very strict electron identification re- 
quirements. Those electrons were paired with all other 
tracks in the event. The invariant mass distribution of 
such pairs is shown by the upper histogram in Fig. 1101 

One can see the characteristic shape of the partially 
reconstructed ir° Dalitz decays and a peak at around 
22 MeV/c 2 corresponding to 7-conversions close to the 
beam pipe. Since the conversion electrons originate at 
the displaced converter vertex, and therefore skip the first 
3.8 cm of the magnetic field, the reconstructed invariant 
mass peak is shifted from zero. Among these pairs a fur- 
ther selection was made to choose those which open up 
in the plane perpendicular to the detector magnetic field. 
This requirement effectively suppresses the combinatorial 
background and pairs coming from the 7r° Dalitz decays, 
but keeps 7 — > e + e~ pairs having small opening angle. 
The middle histogram in Fig. [10] shows that the con- 
version peak significantly dominates the residual Dalitz 
contribution and the combinatorial background. Finally 
we applied the electron identification requirements to the 
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FIG. 8: Detector acceptance as a function of transverse momentum for different particles measured with the PHENIX experi- 
ment. The band in the left panel shows the largest relative systematic uncertainty among all curves. 
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FIG. 9: Efficiencies of the ERT 4x4 and ERT 2x2 triggers for single 7 clusters and electrons as a function of energy (left). 
Trigger efficiencies for different meson decays as determined from Monte Carlo simulation (right). The band shows the largest 
relative systematic uncertainty among all curves. 



second track. The invariant mass distribution of the pairs saturates at <~ 10%. 
where the second track fails to be identified as an electron 
is shown by the filled histogram. The ratio of the lowest 
to the middle histogram under the peak is the electron 
identification loss. It reaches 20% below 0.5 GeV/c and 
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FIG. 10: Invariant mass distribution for e + e~ pairs where 
one track is identified as an electron and the second track is 
any track (dashed line), the same for pairs which open up in 
the plane perpendicular to the detector magnetic field (solid 
line), and among those pairs, the ones in which the second 
track fails the electron identification cut (filled histogram). 



E. Calculation of invariant cross sections 

The invariant cross section for a particle in each pt bin 
was calculated as: 



1 



1 



1 



1 



N(A PT ) 
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(1) 



where N(Apx) is the number of reconstructed particles 
in a given px bin, C is the integrated luminosity sampled 
by the minimum bias trigger, s(j>t) is the acceptance and 
reconstruction efficiency, BR is the branching ratio, and 
£bbc is the minimum bias trigger efficiency for events 
containing mesons, estimated to be 0.79±0.02. The cross 
section sampled by the BBC trigger, a^ t — 23.0±2.2 mb, 
was used to determine the integrated luminosity. For 
the analyses with the minimum bias data sample e(pt) 
corrects for the acceptance and reconstruction efficiency 
while for analyses with the ERT data sample it includes 
the ERT trigger efficiencies as well. A bin shift correction 
was applied to take into account the finite width of the 
Pt bins used in the analyses. The correction is made by 
shifting the data points along the vertical axis according 
to the procedure described in [60} . 

Finally, in the ui —> tt°tt + tt^ and Kg —> ir°ir° analy- 
ses, the cross sections measured with the ERT and with 
the minimum bias triggers were averaged in the overlap- 
ping px region, taking into account the statistical and 
systematic uncertainties. 



estimated by varying the analysis cuts, and by varying 
the energy and momentum scales of the EMCal and DC 
by 1%. The resulting systematic uncertainties for the 
different decay modes of rj, uj, rj', and 4> mesons are 
summarized in Table [I] The uncertainties are categorized 
by types: (A) uncorrelated between pt bins, (B) px cor- 
related, all points move in the same direction but not by 
the same factor, (C) an overall normalization uncertainty 
in which all points move by the same factor independent 
of pt- The type C uncertainty is predominantly due to 
the uncertainty of the minimum bias trigger efficiency in 
p + p collisions, equal to 9.7% [H, 0- The uncertainty of 
the raw yield extraction is estimated as described in sub- 
section IIII Bl It dominates the total uncertainty and is 
split into Type A and Type B contributions. 



G. Neutral meson spectra 

The invariant differential cross sections calculated us- 
ing Eq.[T]are tabulated in Tables IXl and IXII and plotted in 
Fig.rTTJ Different symbols are used to show results for dif- 
ferent decay modes. One can see a very good agreement 
between the particle spectra measured in the different 
decay modes. Results for low pt bins for particles recon- 
structed through decays in the e + e~ mode allow an ac- 
curate measurement of the integrated particle yield. The 
integrated yield at midrapidity for the w is measured to 
be da"/dy = 4.20±0.33 stat ±0.52 s y st mb and for the <t> is 
measured to be da^/dy = 0.432±0.031 stat ±0.051 syst mb. 
The mean transverse momentum for these particles is 
(p T ) = 0.664 ± 0.037 stat ± 0.012 s y st GeV/c and (p T ) = 
0.752 ± 0.032 stat ± 0.014 s y st GeV/c. 



IV. ANALYSIS OF PARTICLE SPECTRA 

In this section we analyze the measured invariant 
transverse momentum spectra for a variety of hadrons in 
p + p collisions at y/s = 200 GeV and search for common 
features. All measurements are quoted as the invariant 
differential cross sections at midrapidity averaged over 
the rapidity interval \y\ < 0.35. 



d 3 a _ inol 



where af^ el = 42 mb. 



1 



1 



d 2 N 



2tt Pt N e . 



dydpj 



(2) 



A. Data samples 



F. Systematic uncertainties 

In addition to the systematic uncertainties described in 
the corresponding analysis sections, uncertainties of the 
ERT trigger efficiency and acceptance corrections were 



The procedures used for the reconstruction of the par- 
ticle transverse momentum spectra are described above 
in section IIIII and in other PHENIX publications listed 
in Table |TTJ Figure [TT] shows (left) the results and (right) 
the results compared with previous PHENIX results. All 



14 



TABLE I: Relative systematic uncertainties (in percent) for different decay modes. Given ranges indicate the variation of the 
systematic uncertainty over the pr range of the measurement. 



Particle 


K°. 




CO 




77/ 


4 




Uncertainty 


Decay 


n°n° 


7r°7r + 7r 


v> 7 


e + e 


rj7V + TT~ 


K+K~ 


e + e 


Type 


Acceptance 


8 


5 


6 


5 


5 


5-7 


5 


B 


EMCal energy resolution 


4-5 


2-5 


2-3 




2-4 






B 


EMCal, DC scale 


4-6 


2-6 


3-17 


2-11 


2-5 


1-5 


2-10 


B 


7T°, 77 selection 


5-10 


3 


3 




3 






B 


ERT trigger efficiency 


2-12 


3-10 


2-7 


1-3 


2-4 




1-2 


B 


Peak extraction MC 


2 


1 


1 


1 


1 


3 


1 


A.B 


Raw yield extraction 


4-19 


5-17 


5-12 


4-15 


6-25 


8-25 


3-11 


A,B 


7-conversion 


6 


3 


5 




3 






C 


e-identification 








10 






9 


B 


Branching ratio 





1 


3 


1.7 


3 


1 


1.3 


C 


MinBias Trigger 


9.7 


9.7 


9.7 


9.7 


9.7 


9.7 


9.7 


c 


Total 


17-29 


13-24 


15-26 


16-24 


14-29 


14-28 


15-18 






2 4 6 8 10 12 

P T (GeV/c) 

FIG. 11: Invariant differential cross-section of neutral mesons 
measured in p + p collisions at i/s =200 GeV in various decay 
modes. The lines are fits to the spectra as described further 
in the text. 



meson spectra used in this paper are not corrected for 
feed down. 

Figure [12] demonstrates a very good agreement be- 
tween the new results and previously published data. 
The results presented in this paper greatly enhance the 
Pt range of the previously measured particles and add 



TABLE II: Data samples used in the analysis of particle spec- 
tra. The [X] and IXII in the "Ref." column refer to Appendix 
Tables IXl and IXB 



Particle 


Mode 


Physics 
Run 


Pt (ttit) range 
GeV/c,GeV/c 2 


Ref. 




7T 


77 


5 


0.5-20 


[3] 


7T + , 7T~ 


ToF 


3 


0.3-2.7 




K+, K~ 


ToF 


3 


0.4-1.9 




K°s 


_0_0 
7T 7T 


5 


2-13.5 




V 


77 


3 


2-12 




V 


77 


6 


2-20 




1 




3 


2.5-8.5 




u 


e + e~ 


5 


0-4 


m 


u 


7r°7r + 7r - 


5 


2-13.5 


x 


Ul 




3 


2.5-10 




u 




5 


2-12 


m 


UJ 


7T°7 


3 


2-7 


[4] 


n' 




5 


3-11.5 


El 


<t> 


e + e~ 


5 


0-4 


IU 


<P 


K + K~ 


5 


1-8 


El 


J/4> 


e + e~ 


5 


0-9 




J/H> 


e + e~ 


6 


0-9 


[48j 


1>> 


e + e~ 


6 


0-7 


M 


p,P 


ToF 


3 


0.6-3.7 





results for particles that have not been previously ana- 
lyzed. 

For each particle we considered all available measure- 
ments of the invariant momentum distributions together 
with their statistical and systematic uncertainties catego- 
rized as Type A, B, and C, as explained in section UlI Fl 

For the analysis of the shape of the transverse momen- 
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FIG. 12: Invariant differential cross sections of different par- 
ticles measured in p + p collisions at y/s —200 GeV in various 
decay modes. The spectra published in this paper are shown 
with closed symbols, previously published results are shown 
with open symbols. The curves are the fit results discussed 
in the text. 



turn distributions the data for all particles of the same 
isospin multiplet were combined into one px spectrum to 
be fitted. All data for positively and negatively charged 
particles measured in the same analysis and in the same 
Pt bins were averaged. All data for neutral particles, 
measured via different decay channels, were added to- 
gether. The notation tt is used to denote a combined 
spectrum of tt° and (n + +ir~)/2, K is used for a combined 
spectra of K s and (K + +K~)/2, p denotes (p+p)/2, and 
so forth. Independent measurements of the same particle 
performed using different data samples or different decay 
modes were also added together but not averaged. For 
data samples where the results were published as dN/ dpT 
a conversion was made using Eq. [2l 



B. Particle spectra fit distributions 

It is widely known from experimental data that, as ex- 
pected from pQCD calculations (e.g. 62]), a pure power 
law shape successfully describes the high px region of 
particle spectra: 



E 



fa 

dp 3 



Ap, 



(3) 



where the shape is determined by the power v and A is 
a normalization constant. However, the power law shape 
is seen to fail in the region below about pr = 3-5 GeV/c 
where the spectra exhibit a more exponential shape. 

The exponential shape of the particle spectra at low 
Pt suggests a thermal interpretation in which the bulk 
of the produced particles are emitted by a system in ther- 
mal equilibrium with a Boltzmann-Gibbs statistical de- 
scription of their spectra: 



E dp^ =Cbe 



(4) 



where C& is a normalization factor and E is the particle 
energy. At midrapidity one can replace E by tot = (Py + 
to 2 ,) 1 / 2 , where mo is the particle rest mass. 

In recent years a variety of publications 0, [2H-I30I l63l | 
have used the Tsallis distribution [24[ to fit particle spec- 
tra. The Tsallis distribution derives from a generalized 
form of the Boltzmann-Gibbs entropy and is written as: 



G q (E) = C q 1 -(!-<?) 



T 



1/(1-9) 



(5) 



q , E, and T have similar meanings as in Eq.|4]and 



where C, 

q is the so-called nonextensivity parameter. For values of 
q =/= 1 the distribution exhibits a power law behavior with 
power n — —1/(1 — q). In order to associate the Tsallis 
distribution with a probability distribution, which de- 
scribes the invariant particle spectra given by Eq. [5] and 
defined over < E < 00, Eq.[5]must satisfy a normaliza- 
tion and energy conservation condition (E) < 00. This 
limits the range of the parameter q to 1 < q < 1^. The 
Tsallis distribution reduces to the Boltzmann-Gibbs dis- 
tribution of Eq. [4] in the limit of q — > 1 . 

To put Eq. [5] into a form appropriate to fit particle 
spectra we replace E by tot = (p\ + to 2 ,) 1 / 2 and use 
the requirement of unit normalization to determine the 
coefficient C q in Eq. [5] to be equal to: 



C 



(2g-3)(«-2) 



T(T + too) -(g-l) (9 -2)to 2 
1 

X (1 - (1 -g)IHll) 1/(1 - 9) ' 



Replacing the parameter q with: 



1 



1-9 



(0) 



(7) 



The resulting formula used in the fitting procedure is 
given by: 



E 



dp 3 



1 da 



(n- l)(n-2) 



27T dy (nT + mo(n — l))(nT + too) 
nT + tot n 



nT + too 



(8) 



where da /dy is the integrated cross section of the particle 
production at midrapidity. 
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In the limit of tuq —> Eq. [8] becomes: 



E 



d 3 a _ 1 da (n - l)(n - 2) 
dp 3 27r dy (nT) 2 



1 



nT/ 



(9) 



This form is very similar to the QCD inspired expression 
suggested by Hagedorn in [22[ written as a function of 
tot instead of pr- 

The condition that the shape of the TOr-spectra of dif- 
ferent particles are the same regardless of their mass, is 
referred to as my-scaling. m-r scaling is known to provide 
a good description of the experimental data at low en- 
ergy, where the spectral shapes are exponential [l8|, [23[ . 
Due to the explicit too mass dependence in Eq. [8] the 
Tsallis distribution docs not satisfy tot scaling, except 
in the case too — > or q — > 1, in which case the limiting 
forms of Eqs. [4] or [9] apply. Therefore the accuracy of fits 
to the Tsallis distribution and the validity of tot scaling 
needs to be quantified with data. 

The power law behavior at high px which appears in 
Eq. [8] is governed by the parameter n. The parameter 
n can be related to the simple power law parameter v 
that occurs in Eq. [3] through the condition that both 
expressions have the same power-like slope at a given 
Pt- From Eqs. [8] and [3] one can write: 



d\n(nT + Tot)" 
d\n(p T ) 



dh\(p T v ) 
dh\(p T ) 



Pt 



vTmi 



(10) 



At high pt (pt 3> ma^T) where one can neglect the 
difference between tot and pt, v and n coincide. In the 
Pt region where most particle spectra are measured, n is 
15-25% larger than v. 

The mean of the Tsallis distribution in the form 
of Eq. |S]is calculated as: 



2nT 



{m T ) 



+ 



(n-2)(n- 1) 



n-3 (raT + m (n- l))(n-3) 
2nT n-2 



n - 



n ■ 



-m . 



(11) 



The approximate relation requires too 3> T. This condi- 
tion is satisfied for all particles, except pions, for which 
T and too are about the same. Similarly, the mean pt 
can be well approximated for all measured particles with 
a linear dependence: 



(Pt) 



2nT 

n — 3 



f(n)m . 



(12) 



The first contribution is identical to that in Eq. [TT] and 
f(n) has only a weak dependence on mo which we neglect 
in Eq. QJ] 

The Tsallis distribution is appealing to use to describe 
particle spectra because it provides a single functional 
form that can reproduce the full spectral shape with just 
two parameters that potentially have an underlying phys- 
ical interpretation. Tsallis distributions have been used 



successfully to describe particle spectra in different colli- 
sion systems and at different energies @, [26l - [30l |63 - (66| . 
Tsallis distributions also describe various physics phe- 
nomena beyond particle production and have been suc- 
cessfully applied in other fields of science, see [HI, l67l - [69| 
and references therein. 

As mentioned above, the Tsallis distribution was de- 
rived as the single particle distribution corresponding to 
a generalization of the Boltzmann-Gibbs entropy through 
the introduction of the non-extensivity parameter q [24| . 
Whereas the Boltzmann-Gibbs distributions are found 
to apply to systems which exhibit an exponential relax- 
ation in time to a stationary state characterized by ex- 
ponentials in energy at thermal equilibrium, the gener- 
alized form is found to apply to systems which exhibit 
power laws in relaxation time and energy. These are sys- 
tems which relax with a non-ergodic occupation of phase 
space as a consequence of the microscopic dynamics of 
the system. Among other examples, this is characteris- 
tic of systems with long range interactions that fall off 
with distance with a power smaller than the dimension- 
ality of the system. It is an interesting question whether 
strongly interacting partonic matter might also exhibit 
power law relaxation. In fact, an analysis of the diffusion 
of a charmed quark in partonic matter produced in par- 
ton cascade calculations found that the parton densities 
were characterized by Tsallis distributions, rather than 
Boltzmann-Gibbs distributions [70l |. 

The physical interpretation of the parameter T in 
Eq. [5l especially in p + p collisions, is not straightfor- 
ward. One can expect that for larger systems, as pro- 
duced in relativistic heavy ion collisions, it reflects the 
kinetic freeze-out temperature (7kf ) at which particle 
scattering ceases to modify the spectral shapes. It is 
shown below that the magnitudes of (T) found in this 
work are close to (Tkfo) extracted in the blast-wave model 
approach @,[7l| applied to p + p data. In high energy ap- 
plications it has been shown 25] that the parameter q of 
the Tsallis distribution of Eq. [5] can be related to the 
amount of temperature fluctuations in the system as: 



Var(jy 

(t) 2 



(13) 



C. Fitting procedure 

In order to obtain a reliable estimate of the fit uncer- 
tainties, the experimental systematic uncertainties must 
be treated properly. The various types of systematic 
uncertainties have been taken into consideration as de- 
scribed here. The pr independent systematic uncertain- 
ties of Type A have been combined in quadrature with 
the statistical errors and the pT-independent systematic 
uncertainties of Type C was reduced by 9.7% due to 
the trigger uncertainty, common to all analyzed particles. 
Residual uncertainties of Type C and of Type B must also 
be considered in the analysis. The Type B uncertainties 
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by definition have an unknown px dependence. In order 
to estimate their effect, the particle spectra were varied 
and fit multiple times. For each fit the y-coordinate in 
each px bin was varied by the same amount according to 
the uncertainty of Type C, and by differing amounts ac- 
cording to the Type B uncertainties, in a manner similar 
to that explained in [72j . 

Variations of the y-coordinates were made indepen- 
dently for each fit with the amount of variation chosen 
randomly according to the p^-dependcnt uncertainties 
for each particle and each sample. For the particle spec- 
tra consisting of multiple samples results of each fit to 
the entire spectrum were weighted with the probability 
of the fit estimated from the x 2 criteria. Such weighting 
emphasizes variations in which individual samples fluctu- 
ate toward each other rather than away from each other, 
which corresponds to the assumption that the different 
samples represent measurements of the same true mo- 
mentum distribution. 

As a result of the multiple fits, weighted distributions 
of the fit parameters were obtained. The mean of the 
distribution was taken as the parameter value, the RMS 
width of the distribution was taken as the systematic un- 
certainty, and the statistical uncertainty was taken from 
the fit to the unmodified data. The number of fits was 
chosen such that the mean and the RMS did not change 
with increasing number of trials. 

D. Fit results 

The fits of Eq. [8] to the data are shown in Fig. [T3l with 
dotted lines. The results are given in Table HTT1 

TABLE III: Parameters of the Tsallis fit with Eq.|B]with all 
parameters free to vary. The uncertainties are statistical and 
systematic. Cross sections are in fib for J/ip and ip' and in 
mb for all other particles. 

da/dy (mb, fib) T (MeV) n = -1/(1 - q) 

7T 43.5 ±2.0 ±1.9 112.7 ±2.9 ±1.1 9.57 ± 0.11 ± 0.03 

K 4.0 ±0.1 ±0.5 132.7 ±3.8 ±7.2 10.04 ± 0.16 ± 0.27 
77 5.1 ±1.1 ±3.9 119 ±10 ±30 9.68 ± 0.18 ± 0.49 

w 4.3 ±0.3 ±0.4 109.7 ±6.9 ±6.7 9.78 ± 0.24 ± 0.18 

T}' 0.80 ±1.5 ±0.7 141 ±107 ±61 10.5 ±2.2 ±1.2 
4> 0.41 ±0.02 ±0.03 139 ±16 ±15 10.82 ± 0.71 ± 0.56 
J/ip 0.73 ± 0.01 ± 0.05 149 ± 56 ± 82 12.3 ± 1.6 ± 2.9 

Tp' 0.13 ±0.03 ±0.02 164±10 3 ±10 2 14 ±12 ±6 
p 1.63 ±0.05 ±0.11 107 ±13 ±12 12.2 ±1.0 ±0.7 



The fit parameters n and T are strongly correlated. In 
some cases, the parameter T can change by more than 
a factor of two, and still produce a good fit. Therefore 
additional information is needed to constrain the values 
of n and T. For that purpose one can use a power law 
given by Eq. [3] fitted to the same data. As discussed 
above the parameters n and v are related to each other 
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FIG. 13: The pr spectra of various hadrons measured by 
PHENIX fitted to the power law (dashed lines) and Tsallis fit 
(solid lines). See text for more details. 

through Eq. QUI However, it is found that the results of 
the power law fit depend on the fit range, but become 
stable when the fit range begins above pt ~ 3.5 GeV/c 
for most particles, or above pr ~ 5.5 GeV/c for heavier 
particles such as the J / ip. The resulting power law fits are 
shown in Fig. [T3] as dashed lines that have been plotted 
down to pt =0.5 GeV/c. Spectra without sufficient data 
above the fit range lower limit were not fitted. The results 
are given in Table IIVI 

TABLE IV: Parameters of the power law fit with Eq. |3] The 
uncertainties are statistical and systematic. Units of A are 
mb(GeV/c)"+ 2 . 



v A 



7T 


8.174 ±0.035 ±0.049 


16.4 ± 1.1 ± 1.6 


K 


8.24 ±0.08 ±0.11 


8.8 ± 0.9 ± 1.6 


V 


8.169 ±0.037 ±0.054 


7.64 ± 0.46 ± 0.83 




7.986 ± 0.083 ± 0.080 


9.5 ± 1.3 ± 1.4 


rf 


8.12 ±0.21 ±0.11 


3.6 ± 1.2 ±0.8 




8.20 ±0.36 ±0.15 


2.8 ± 1.5 ±0.7 




7.0 ±1.2 ±0.4 


0.03 ±0.03 ±0.02 



The parameters v of the power law fits and the pa- 
rameters n and T of the Tsallis fits are shown in Fig. [TJ] 
as a function of the particle mass. The parameters have 
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been fit to a linear function to establish if there is a mass 
dependence. The fits are shown in Fig. Q3] as solid lines 
with the uncertainties indicated by dashed lines. From 
Fig. [TH it is evident that the parameters are consistent 
with no significant mass dependence. Therefore the pa- 
rameters have also been fit with a constant value. The 
results for the linear and constant fits are summarized in 
Table. El 
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FIG. 14: Particle mass dependence of the fit parameters. 
Power law parameters v and n are plotted in the upper panel. 
Vertical bars denote the combined statistical and systematic 
uncertainties. The solid lines are linear fits. The dashed lines 
denote the fit uncertainty within which the linear fit can be 
inclined. The lower panel shows the same for the fit param- 
eter T. The proton measurement (open circle) is not used in 
the fits. 



TABLE V: Constant and linear fits to the power law and 
Tsallis fit parameters. The last column (Prob.) gives the 
probability estimated by the y 2 /n.d./. of the fit." 





Fit 


Prob 


v 


8.154 ±0.039 


0.75 


v 


(8.22 ± 0.07) - (0.15 ± 0.14)m [GeV/c 2 ] 


0.79 


n 


9.656 ± 0.097 


0.69 


n 


(9.48 ± 0.14) + (0.66 ± 0.39)m [GeV/c 2 ] 


0.94 


T (MeV) 


115.3 ±2.8 


0.43 


T (MeV) 


(111.5 ± 4.0) + (15 ± 12)m [GeV/c 2 ] 


0.51 



The fitted linear coefficients are consistent with zero 
within less than two standard deviation of the fit ac- 
curacy for all three parameters. At the same time the 
parameter v is more accurately defined compared to the 
Tsallis fit parameter n. We can invoke Eq.[TU]to constrain 
the Tsallis fit using the parameter v. This requires to es- 
timate the effective px which appears in Eq. 1101 Using 
the mass independent terms of the fits listed in Table |V] 
the effective pr is about 7 GeV/c. 



This value is large enough to allow to neglect the dif- 
ference between mx and pt in Eq. [TU] for all particles, 
except the J/tp and ip'. These two particles do not con- 
strain the mass dependence of the Tsallis fit parameters 
due to their large fit uncertainties, as shown in Fig. I14I 

Under the assumption that the parameter v is the same 
for all particles, the mass dependence of the parameters 
n and T must either be present or absent together. This 
can be checked by fixing the parameter n to a constant 
value of n — 9.656 (from Table fV|) and fitting the data 
again. The mass dependent coefficient for the parameter 
T that results in this case is somewhat different from zero 
compared to uncertainties. This is a clear contradiction 
to Eq. II Pl under the assumption of constant v, and there- 
fore indicates that the parameters n and T have a mass 
dependence. However, this conclusion is at the limit of 
the accuracy of the currently available data. 

For further analysis the parameter n was fixed to have 
a linear dependence n — 9.48 + 0.66too [GeV/c 2 ] (from 
Table [Vj and the particle spectra were fit again. The 
results are given in Table I VII and the fit to the mass 
dependence of T is given in Table IVHI 

TABLE VI: Parameters of the Tsallis fit with Eq.[S]with pa- 
rameter n constrained to a fixed linear dependence on mass 
(for mesons). The uncertainties for da/dy and T are statis- 
tical and systematic, and are only systematic for n. Cross 
sections are in fib for J /if) and tp' , and in mb for all other 
particles. 





da/dy (mb, /xb) 


T (MeV) 


n= -l/(l-g) 




42.8 ±3.1 ±2.7 


112.6 ±2.1 ±2.8 


9.57 ±0.10 


K 


4.23 ±0.09 ±0.53 


125.4 ±0.9 ±5.3 


9.81 ±0.13 


n 


3.86 ±0.30 ±0.71 


124 ± 2 ± 12 


9.84 ±0.14 


w 


4.26 ±0.23 ±0.33 


115.5 ±2.1 ±6.8 


10.00 ±0.22 


rf 


0.63 ±0.27 ±0.21 


123 ± 17 ± 18 


10.12 ±0.28 


4> 


0.427 ±0.019 ±0.023 


123.4 ±3.0 ±8.3 


10.16 ±0.31 


J/tp 0.760 ± 0.014 ± 0.048 


148 ± 8 ± 35 


11.5 ±1.1 




0.132 ± 0.029 ±0.020 


147 ± 127 ± 54 


11.9 ± 1.3 


V 


1.775 ± 0.044 ± 0.066 


58.8 ± 1.8 ±6.1 


9.20 ±0.28 



Comparison of the results listed in Tables|TTT]and|Vl]re- 
veals that the parameters of the fit did not change signif- 
icantly within uncertainties, even for the rj and rj' mesons 
which are not measured at low px- In addition, with the 
n parameter constrained the uncertainty on the parame- 
ter T is reduced. 

Since there is not yet a published PHENIX measure- 
ment of protons at high-p^ the parameter v cannot be 
determined for the case of protons. Results published 
in Q suggest that the slope of the proton spectra at high 
Pt is the same as that for mesons. Using this assump- 
tion allows to extract the parameter T for protons with 
the result listed in Table IVII The value of T for protons 
differs from the values extracted for mesons. 

Using the linear dependence of the T parameter T = 
112.6 + 11.8m [GeV/c 2 ] extracted from the fits to the 
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TABLE VII: Constant and linear fits to the Tsallis parameter 
T of mesons with fixed parameter n. The last column (Prob.) 
gives the probability estimated by the \ 2 /n.d.f. of the fit. 





Fit 


Prob. 


T (MeV) 


117.4 ±2.5 


0.64 


T (MeV) 


(112.6 ± 3.8) + (11.8 ± 7.0)m [GeV/c 2 ] 


0.83 



Tsallis distribution with fixed linear dependence of the 
n parameter (from Table IVI|) the spectra can be fit once 
again to obtain an improved normalization parameter. 
The resulting fits are shown in Fig. [13] as the solid lines, 
and the results of the fit are given in Table I Villi 



TABLE VIII: Parameters of the Tsallis fit with Eq. [H] with 
parameters n and T constrained to have a fixed linear depen- 
dence on mass (for mesons). The uncertainties for da/dy are 
statistical and systematic, and are only systematic for T and 
n. Cross sections are in fib for J/tj) and ip' , and in mb for all 
other particles. 





da/dy (mb, fib) 


T (MeV) 


n = -l/(l-5) 




40.5 ±0.3 ±5.8 


114.2 ±4.0 


9.57 ±0.10 


K 


4.71 ±0.06 ±0.48 


118.4 ±5.2 


9.81 ±0.13 


n 


4.46 ±0.05 ±0.97 


119.0 ±5.4 


9.84 ±0.14 




3.64 ±0.07 ±0.77 


121.8 ±6.7 


10.00 ±0.22 


rf 


0.62 ±0.04 ±0.16 


123.8 ±7.7 


10.11 ±0.28 


<\> 


0.421 ± 0.009 ± 0.054 


124.5 ±8.1 


10.15 ±0.31 


J/1> 


0.761 ±0.013 ±0.060 


149 ± 22 


11.5 ± 1.1 


i>> 


0.133 ±0.024 ±0.019 


156 ± 26 


11.9 ± 1.3 


V 


1.76 ±0.03 ±0.16 


58.8 ±6.4 


9.20 ±0.28 



The parameters n and T, and their errors, are fixed to 
the values obtained from the fitted linear dependence of 
the parameters on particle mass, obtained from the fits of 
TableslVland lVIll The systematic error on the integrated 
yields reflects the variation of the n and T parameters 
within the errors. It also includes the uncertainty from 
the variation of the spectral shapes within errors of Types 
B and C, as explained above. 

The fits accurately describe the data. To demonstrate 
the quality of the fits the data points have been divided 
by the fit value and the ratios plotted in Fig. [T5l 

Grey error bars show the combined systematic uncer- 
tainty of Type B and C, with the Type B uncertainties 
dominating. The dashed lines show the fit uncertainty 
corridor. The RMS of the vertical spread of all points 
plotted in Fig. [T5]is 0.17. If each point is normalized to 
the combined statistical and systematic error of the data 
point, the RMS of the same distribution is much larger 
with a value of 0.88, which indicates that the agreement 
between data and fit is well within errors. 



V. DISCUSSION 

A. Tsallis fit parameters 

The analysis of section IIVI demonstrated the ability 
of the Tsallis distribution functional form to fit the full 
transverse momenta spectra for all different species pro- 
duced in p + p collisions at y/s = 200 GeV with only two 
parameters, n = —1/(1 — q) and T. Furthermore, the 
values of the two parameters extracted from the fits are 
approximately the same for all measured mesons. 

On the other hand, the observation that the pure 
power law fit of Eq. [3] to the spectra in the region of 
Pt >3.5 GeV/c yields the same power v = 8.154 ± 0.039 
for all particles with higher accuracy than the Tsallis 
fit, indicates that a weak mass dependence of the Tsal- 
lis parameters is to be expected. Assuming a weak 
mass dependence one gets T = 112.6 ± 3.8 ± (11.8 ± 
7.0)m [GeV/c 2 ] MeV and n = 9.48 ± 0.14 ± (0.66 ± 
0.39)m [GeV/c 2 ] that improves the description of the 
meson spectra with the Tsallis distribution. The pa- 
rameters are listed in Tables IIVI and I Villi and plotted 
in FigQH 

The ratios of the data points to the Tsallis parame- 
terization using the global fit parameters n and T for all 
particles were shown in Fig. 1151 The figure represents 
nine different particles species measured over the range 
< pt (GeV/c) < 20 using six independent data samples 
and ten different analysis techniques. The parameteriza- 
tion is in good agreement with the experimental data. 
The average deviation of the points from one in all pan- 
els of the figure is 88% of the combined uncertainty of 
the data and the fit. 

The Tsallis distribution fit for the proton measurement 
yields a parameter T — 58.8±6.4 MeV significantly lower 
than that for the mesons. Since the published PHENIX 
results for protons have limited pr range this result was 
checked and confirmed using STAR measurements for 
protons and heavier baryons 04H1- This result indi- 
cates significantly different Tsallis fit parameters between 
mesons and baryons. 

The similarity of the measured parameters T and n for 
all studied mesons suggests a similar production mecha- 
nisms in p ± p collisions at s/s — 200 GeV. At the same 
time, the mechanism of baryon production must have dif- 
ferent features. The interpretation of the T parameter of 
the Tsallis fits is not straightforward. If interpreted as 
a temperature the values obtained are seen to be similar 
to average freeze-out temperatures (Tkfo) extracted in the 
blast- wave model approach [8j, l71| applied to p ± p data. 
As mentioned above, the parameter n can be related to 
temperature fluctuations as y/Var(l/T)/(l/T) = 1/n in 
a thermal interpretation. Following this interpretation 
one can estimate the fluctuations of the inverse slope pa- 
rameter 1/T to be of order of 0.3. 
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FIG. 15: Data to fit ratio for different particies used in the anaiysis. The systematic uncertainties are the combined uncertainties 
of Type B and Type C, excluding the common 9.7% trigger uncertainty. 



B. vtiT scaling 

As discussed in section IIVB1 mr scaling can not be 
an exact scaling when particle spectra follow the Tsallis 
distribution with g ^ 1. However, rriT scaling might 
be found to be approximately true. The validity of mr 
scaling can be studied quantitatively with the assistance 
of Eq. IH1 which gives the Tsallis distribution in the limit 
too — > with a form explicitly satisfying scaling. 

Figure [16] shows the spectra for all particles plotted 
as a function of tot and normalized at one single point 
on the X-axis. All normalized spectra are then fit simul- 
taneously with Eq. [5] using fixed parameters taken from 



Tables EIVnin=9. 656 and T = 115.3 MeV for mesons, 
and T = 58.8 MeV for baryons. 

The difference in the spectral shapes between mesons 
and baryons is apparent from the figure. It is due to the 
large difference in (T) between these particle groups. At 
the same time the spectra of both mesons and baryons 
separately are well described by the mj- scaling assump- 
tion. 

To quantify this statement we restricted the analysis to 
the PHENIX meson measurements only. After optimiza- 
tion of the normalization point for the different particles 
the RMS of the data to fit ratio for all points shown in 
Fig [TB] has a value of 0.25. This is to be compared to 
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1 m T (GeV/c 2 ) 10 

FIG. 16: Particle spectra plotted vs. rnr and arbitrarily nor- 
malized at pt =10 GeV/c. Open symbols are mesons and 
full symbols are baryons measured by PHENIX (circles) and 
STAR (stars). The lines are the Hagedorn fits by Eq. [5] to 
mesons (solid) and baryons (dashed) with parameters n and 
T fixed to average values. Error bars are statistical and point- 
by-point systematic only. 

the analogous result of Fig. [15] for the Tsallis fit in px 
coordinate which gave an RMS of 0.17. This small in- 
crease supports the conclusion that at \/s =200 GeV all 
meson spectra have very similar shape when plotted as a 
function of my, and thus obey raj- scaling. 



C. Integrated yields and (pr) 

Using the Tsallis functional form and Tables fVIVIII one 
can derive information about (ttit) and (jot), based only 
on the particle mass and the baryon number. Deter- 
mination of the integrated da/dy requires experimental 
measurement of the particle production cross section in 
at least in a limited pr range. 

The results presented below were obtained indepen- 
dently for each particle species without averaging within 
the same isospin multiplet, unless such averaging was 
done by the experiment. Different measurements of the 
same particle were combined together. Published data 
from the STAR experiment and the references listed in 
Table iFXl were also analyzed. To compare PHENIX and 
STAR results, the spectra and the integrated yields pub- 
lished by STAR, in units of dN/dy, were multiplied by 



30 mb, which is the value of the STAR minimum bias 
cross section in p + p collisions at y/s — 200 GeV, that 
includes the non-single diffractive part of p + p interac- 
tions (cf. 0). 

The particle spectra published by the STAR exper- 
iment were fit to the Tsallis functional form given by 
Eq. U with the parameters n = 9.48 + 0.66ra [GeV/c 2 ] 
and T = 112.6 + 11.8m [GeV/c 2 ] taken from the global 
fit to the PHENIX data. The same parameters deter- 
mined independently for the STAR data give consistent 
results for mesons. For baryons the STAR data showed 
a dependence of the parameter T on the mass of the 
particle, however the fit uncertainties were too large to 
make a definite statement. The value of T averaged over 
all baryon measurements made by STAR agrees with the 
PHENIX result for the proton measurement. Calculation 
of da/dy for p and p measured by PHENIX was not done 
because the spectra are feed-down corrected and the ex- 
trapolation to low px requires additional evaluation of 
the systematic uncertainties. 

Figure [17] shows a comparison of the experimentally 
measured integrated spectral characteristics to the re- 
sults obtained using the Tsallis fits. The ratio of the 
measured characteristic width to the width calculated 
from the Tsallis fit is shown in the upper panel. For 
most particles the width is taken to be (pt), but for the 
J/ip and ip' the comparison is done for (p T ) because this 
is the parameter published in the corresponding articles. 
Statistical and systematic uncertainties of the published 
results are shown at each data point and the uncertain- 
ties of the Tsallis fit values are shown by the band around 

y = i- 

For all mesons the agreement between the published 
values and the values from the Tsallis fit analysis is con- 
sistent with the published uncertainties. This demon- 
strates the accuracy to which the Tsallis functional form 
describes the experimental spectral shapes. 

Eq- HU suggests that the mass dependence of the (j?t) 
should be approximately linear. A fit to the average mean 
momentum of all mesons extracted from the Tsallis dis- 
tribution fits as a function of their mass gives (jpr) — 
(0.319 ± 0.007) [GeV/c] + (0.491 ± 0.009)m . A fit to the 
published data directly gives a similar consistent result of 
(p T ) = (0.284 ± 0.015)[GeV/c] + (0.506 ± 0.033)m . For 
baryons the agreement with the linear fit is reasonable 
based on the data published by the STAR experiment. 

In the original work of R. Hagedorn a nearly linear 
dependence of the (pr) was derived based on the assump- 
tion of Boltzmann-Gibbs statistics to describe the parti- 
cle spectra at low pt- The difference between mesons and 
baryons would follow from the bosonic and fermionic na- 
ture of these particles. However, quantitatively the val- 
ues of particle (jpt) &ud the magnitude of the meson-to- 
baryon difference are not the same as would follow from 
the mechanisms discussed in [22j ■ 

The lower panel of Fig. [T7] shows the ratio of the in- 
tegrated yields published by the experiment to the inte- 
grated yields extracted from the Tsallis function fits. The 
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TABLE IX: Cross sections in mb and (pr) in GeV/c of different particles in p + p collisions at y/s =200 GeV. PHENIX and 
STAR columns show the values obtained by fits to the experimental spectra with the Tsallis functional form as described in the 
text. One should state explicitly that these values do not supersede values given in the "Published" column by the experiments, 
in their publications listed in the last column, or elsewhere. An additional 9.7% systematic uncertainty should be added to 
all da/dy values listed in the column "PHENIX" and 12% to the values in the column "'STAR' to account for the trigger 
uncertainties. Values in the column "Published" are also given without these systematic uncertainties. The column "S.M." is 
the prediction of the statistical model discussed in the text. The characteristic widths of the particle spectra are (pr) for all 
species except for J/ip and ip' for which the values given in the Table are (Pr)- For if)' the integration is done in the pr region 
below 5 GeV/c. All errors are the combined statistical and systematic uncertainties. 

da/dy (mb) (p T ) (GeV/c), ( P 2 T ) (GeV 2 /c 2 ) 



Particle 


PHENIX 


STAR 


Published 


S.M. 


Fit 


Published 


Ref. 


7T° 


41.4 ±5.8 






46.9 


0.377 ±0.012 






7T+ 


39.4 ±7.3 


43.8 ±3.3 


43.2 ±3.3 


42.1 


0.379 ±0.012 


0.348 ±0.018 


[8] 


TT~ 


38.6 ±7.2 


43.2 ±3.3 


42.6 ±3.3 


41.5 


0.379 ±0.012 


0.348 ±0.018 


f8] 


K+ 


4.57 ±0.61 


4.72 ±0.39 


4.50 ± 0.39 


4.57 


0.567 ±0.017 


0.517 ±0.030 


f8] 


K~ 


4.20 ±0.51 


4.61 ±0.18 


4.35 ± 0.39 


4.38 


0.567 ±0.017 


0.517 ±0.030 


[8] 


K% 


5.28 ±0.53 


4.26 ±0.15 


4.02 ±0.34 


4.40 


0.569 ±0.017 


0.605 ± 0.025 


f9] 


n 


4.47 ±0.96 






4.93 


0.595 ±0.018 






p 




6.55 ±0.37 


7.8 ± 1.2 


5.58 


0.714 ±0.019 


0.616 ±0.062 


[73] 


id 


3.65 ±0.77 




4.20 ± 0.47 


5.03 


0.718 ±0.022 


0.664 ± 0.039 


this work 


v' 


0.62 ±0.17 






0.365 


0.808 ± 0.026 






(K*+ + K*-)/2 




1.46 ±0.10 




1.57 


0.774 ± 0.022 






(K*° + K*°)/2 




1.525 ±0.091 


1.52 ±0.19 


1.55 


0.776 ± 0.022 


0.81 ±0.14 


[63] 


4> 


0.421 ±0.055 




0.432 ± 0.035 


0.339 


0.839 ± 0.027 


0.752 ± 0.043 


this work 


<t> 




0.525 ±0.018 


0.540 ± 0.086 


0.339 


0.839 ± 0.025 


0.820 ±0.051 


[74] 


J/%1) (xlO 3 ) 


0.759 ± 0.053 




0.746 ± 0.089 




4.464 ± 0.606 


4.60 ±0.19 


® 


ip' (xlO 3 ) 


0.133 ±0.031 




0.126 ±0.034 




4.807 ± 0.443 


4.7 ± 1.3 


[61] 


P 




4.06 ± 0.23 


4.14 ±0.30 


4.47 


0.648 ±0.019 


0.661 ±0.022 


M 


V 




3.28 ±0.23 


3.39 ± 0.36 


3.59 


0.648 ±0.019 


0.661 ±0.022 


f8] 


A 




1.33 ±0.13 


1.31 ±0.12 


1.30 


0.742 ± 0.023 


0.775 ± 0.040 


[9] 


A 




1.20 ±0.12 


1.19 ±0.11 


1.11 


0.742 ± 0.023 


0.763 ± 0.040 


f9] 






0.094 ± 0.020 


0.078 ± 0.028 


0.092 


0.850 ± 0.030 


0.924 ± 0.054 


f9] 






0.091 ±0.019 


0.087 ± 0.031 


0.082 


0.850 ± 0.030 


0.881 ±0.051 


® 


£*+ 




0.358 ± 0.026 


0.321 ±0.044 


0.308 


0.882 ± 0.032 


1.020 ± 0.073 


[10] 


E*+ 




0.310 ±0.025 


0.267 ± 0.038 


0.260 


0.882 ± 0.032 


1.010 ±0.061 


[10] 


A* ± A* 




0.127 ±0.013 


0.104 ±0.017 


0.168 


0.955 ± 0.038 


1.08 ±0.10 


[10] 


or + Q + (xio 3 ) 




11.5 ±4.6 


10.2 ±5.7 


17.1 


1.035 ± 0.046 


1.08 ±0.30 


[9] 



common uncertainties on all integrated yields of 9.7% for 
PHENIX and 12% for STAR are not included. Most of 
the ratios equal 1 within uncertainties. From Fig. 1171 and 
Table IIXI one may conclude that the constrained Tsallis 
fit reproduces the measured integrated cross section with 
high accuracy for all identified particles in p + p collisions 
at y/s — 200 GeV. This gives justification to use the con- 
strained Tsallis fit results to obtain da/dy for particles 
which have only been measured in a limited px range, 
such as 7T°, 77, and 7/ mesons. The resulting da/dy for 
such particles are also given in Table llXl 

It should be noted explicitly that the da/dy and (pr) 
values given in Table iLXl determined using the Tsallis pa- 
rameterization do not supersede, or presume to be more 
accurate than the corresponding values published by the 



experiments in the original papers. They are given to val- 
idate the method. In those cases where no values have 
been published the Tsallis fit result values in the table 
represent a best attempt to obtain the cross section or 
(pr) based on the validity of the Tsallis fit distribution. 



D. Statistical model calculation 

Figure [T51 shows the ratio of the constrained Tsallis fit 
results for the integrated particle yields to the predicted 
yields from a statistical model (SM) calculation [75j . The 
data-to-model ratio for PHENIX data is shown in the up- 
per panel, and for STAR data in the lower panel. The 
statistical model calculation parameters were chosen to 
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FIG. 17: Comparison of the integrated parameters of the par- 
ticle spectra. The upper panel shows the ratio of the pub- 
lished result divided by the result of constrained Tsallis fit 
for (pt) {{Pt) f° r the J/i/> an d the same for ip' integrated in 
the pt range below 5 GeV/c). The lower panel shows the 
ratio for da/dy. The statistical and systematic uncertainties 
shown at each point are from the published data only. The 
band around y — 1 shows the uncertainty of the values ex- 
tracted from the Tsallis function. The trigger efficiency scale 
uncertainty of 9.7% for PHENIX and 12% for STAR is not 
plotted. There are no points outside the plot boundaries. 
Vertical dashed line separates mesons and baryons. 



FIG. 18: Ratio of integrated yields predicted by the statistical 
model ASM t° those of the constrained Tsallis fits for various 
particles. Results for fits to PHENIX data are shown in the 
upper panel and STAR data in the lower panel. The band 
reflects the uncertainty of the Tsallis fit results and includes 
the trigger uncertainty of 9.7% for PHENIX. The lower panel 
has smaller uncertainties because the model prediction was 
based on published STAR data given in Table IIXI There are 
no points outside the plot boundaries. Vertical dashed line 
separates mesons and baryons. 



VI. SUMMARY 



reproduce the integrat ed y ields published by the STAR 
experiment l63l . 1 73l |74| ]. which may explain the 

larger discrepancies in the comparison to the PHENIX 
results. 

Although statistical models are not commonly used 
to describe p + p data the agreement of the statistical 
model calculation with the STAR results was found to 
be accurate for most particles except for the p, (f>, and 
A* Leaving aside baryons, for which the calcula- 

tions of the da/dy requires additional assumptions, as 
explained above, the Tsallis fit also has difficulty to re- 
produce the result for the p meson as shown in Fig. 1171 
This can be explained by the large systematic uncertainty 
of the published value [73j . 

For the PHENIX data the SM calculations agree with 
the production rates for most mesons because the Tsallis 
fit results of the PHENIX and STAR data agree. The 
production rates of 7r°, 77, ui, rj 1 , and 4> were not measured 
by STAR and so were not used in the determination of 
the SM parameters. Among them the predicted yields 
of 7T° and r\ mesons are in very good agreement with the 
PHENIX data. The prediction of the SM for the w, rf, 
and <p yields are less accurate with ratios just outside of 
errors. 



A systematic study of neutral meson production in 
p + p collisions at y/s —200 GeV has been performed by 
the PHENIX experiment at RHIC with results presented 
in this paper. New measurements by PHENIX of Kg, uj, 
</>, and 77' meson production have been presented. 

The measurement of the K s invariant differential cross 
section via the 7r°7r° decay channel in the momentum 
range 2 < px (GeV/c) < 13.5 extends previously pub- 
lished K ± measurements [l|. 

We present the first measurement of the (f> invariant 
differential cross section in the K + K~ decay mode us- 
ing several different techniques. The combined spectrum 
reaches the upper px limit of 8 GeV/c. 

This work also presents the first measurement of the 
invariant differential cross section of rj' production mea- 
sured via the t]tt + tt~ decay mode with results that cover 
the range 3 < px (GeV/c) < 11. 

Measurements of uj meson production in non-leptonic 
decay channels extends the »t coverage of the previous 
PHENIX uj measurement Q, obtained with a smaller 
data sample, to 13.5 GeV/c. 

First measurements of the uj and <f> in the e + e~ decay 
channel extend the px coverage for these two particles 
down to zero momentum and allow a direct calculation 
of the integrated yields and mean transverse momenta 
with results: da^/dy = 4.20 ± 0.33 stat ± 0.52 s y st mb and 
da^/dy = 0.432 ± 0.031 stat ± 0.051 syst mb; and (p T ) = 
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0.664 ± 0.037 stat ± 0.012 s y st GeV/c and (p T ) = 0.752 ± 
0.032 stat ± 0.014 s y st GeV/c. 

All measured results were found to be consistent be- 
tween the different decay modes and analysis techniques, 
as well as with previously published data. The results 
are shown in Figs. QT] and [12J and the measured cross 
sections are tabulated in the Appendix. 

The invariant differential cross sections for all mea- 
sured hadrons produced in p + p collisions at ^Js 
=200 GeV presented in this work as well as in previ- 
ous PHENIX publications, were shown to be described 
well over the entire momentum range by the Tsallis dis- 
tribution functional form with only two parameters, T 
and n, characterizing the low- and high-pT regions, re- 
spectively. Furthermore, the values of the two parame- 
ters extracted from the fits are approximately the same 
for all measured mesons with a weak mass dependence: 
T = 112.6 ± 3.8 + (11.8 ± 7.0)m [GeV/c 2 ] MeV; and 
n = 9.48 ± 0.14 + (0.66 ± 0.39)m [GeV/c 2 ]. 

The meson spectral shapes have very similar forms 
when plotted as a function of uit and hence follow m-r- 
scaling well at ^fs =200 GeV. On the other hand, the 
proton spectra are described with a significantly lower 
parameter value of T = 58.8 ± 6.4 MeV and do not fol- 
low the mT-scaling form observed for mesons. 

The ability to successfully describe all particle spectra 
in p + p collisions at ^/s = 200 GeV with a common func- 
tional form allows one to calculate the invariant differen- 
tial cross section for any particle. This allows the abso- 
lute integrated yield to be derived from any experimental 
measurement of the hadron spectrum, even with limited 
Pt range. The values of da/dy and (pt) are tabulated in 
Table IIXI for hadrons measured by PHENIX, as well as 
those measured by the STAR experiment using the set of 
values of Tables [V] and I VIII For all measured mesons the 
average transverse momentum of the particle depends lin- 
early on the mass uiq and can be parameterized with the 
relation (p T ) = (0.319±0.007)[GeV/c]+(0.491±0.009)m 

The predictions of statistical model calculations based 
on data published by the STAR experiment [io| were 
shown to be in good agreement with the integrated yields 
calculated from the Tsallis distribution fits for most par- 
ticles. Some deviations are seen for the u>, the rj' , and 
the <f> mesons. 
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1 rl^ 

TABLE X: The invariant differential cross section ^— -±r- of cj meson production measured in the indicated decay channel. 
Notations are: V is the differential cross section, A, B, and C are the three types of errors described in the text. 



meson 


decay 


Pt 


V 


A 


B 






channel 


KjjC V / C 




mb/(GcV/c) 2 








0.125 


3.5 


1.1 


0.4 


0.2 






0.375 


1.76 


0.33 


0.19 


8. x 10 2 






0.625 


1 .12 


0.12 


0.12 


5. x 10 2 






0.875 


0.425 


5.6 x 10~ 2 


4.6 x 10~ 2 


1.9 x 10 2 


UJ 


e + e 


1.125 


0.213 


2.8 x 10~ 2 


2.3 x 10~ 2 


2. x 10 3 






1.375 


9.U X 1U 


1.4 x 10~ 2 


1. x 10~ 2 


4. x 10 3 






1.75 


z. ( o X lu 


4.6 x 10~ 3 


2.9 x 10~ 3 


i i \s in — 3 
l.z X lu 






2.5 


O.o4 X 1U 


6.3 x 10~ 4 


3.5 x 10~ 4 


1.5 X 10 4 






3.5 


3.6 x 10 4 


1.2 x 10" 4 


4. x 10" 5 


2. x 10 5 






2.25 


8.0 x 10 ~~ 3 


1.4 x 10~ 3 


8. x 10~ 4 


2.6 x 10 -4 






2.75 


2.50 X 10 -3 


2.5 x 10~ 4 


2.4 x 10~ 4 


8. x 10~ 5 






3.25 


7.89 x 10~ 4 


5.8 x 10~ 5 


7.4 x 10~ 5 


2.6 x 10~ 5 






3.75 


2.56 X 10~ 4 


1.7 x 10~ 5 


2.3 x 10 -5 


8. x 10~ 6 






4.25 


n a i •. , in — 5 
9.41 X 10 


5.9 x 10~ 6 


8.0 x 10~ 6 


3.1 x 10 6 






4.75 


o an ~.s in — 5 

o.oy X 1U 


2.5 x 10~ 6 


3.5 x 10~ 6 


i o w in — 6 
l.z X 1U 






5.25 


i \ , in — 5 

1.68 X ID 


1.3 x 10~ 6 


1.5 x 10~ 6 


5. x 10 -7 






5.75 


i .0 i X 1U 


7.1 x 10~ 7 


7.1 x 10~ 7 


2.5 x 10 7 






6.25 


o on w in — 6 

3.89 X 1U 


4.1 x 10~ 7 


3.7 x 10~ 7 


1.3 x 10 7 






6.75 


2.13 X 1U 


2.8 x 10~ 7 


2.2 x 10~ 7 


7. x 10 8 


UJ 


_o_+ - 

7T 7T 7T 


7.25 


i A w in — 6 
1.40 X 1U 


2.2 x 10~ 7 


1.6 x 10~ 7 


5. x 10 8 






7.75 


8.5 x 10~ 7 


1.6 x 10~ 7 


1.0 x 10~ 7 


3. x 10~ 8 






8.25 


4.03 x 10~ 7 


9.8 x 10~ s 


4.6 x 10~ s 


1.3 x 10" 8 






8.75 


2.93 x 10~ 7 


7.0 x 10~ s 


3.6 x 10~ s 


1.0 x 10" 8 






9.25 


2.48 x 10~ 7 


6.2 x 10~ s 


3.5 x 10~ s 


8. x 10" 9 






9.75 


1.49 x 10~ 7 


4.0 x 10~ s 


2.2 x 10~ s 


5. x 10" 9 






10.25 


1.09 x 10~ 7 


3.0 x 10~ s 


1.7 x 10~ s 


4. x 10" 9 






10.75 


6.0 x 10~ 8 


1.9 x 10~ s 


1. x 10~ s 


2. x 10" 9 






11.25 


5.1 x 10~ 8 


1.6 x 10~ s 


8. x 10~ 9 


2. x 10" 9 






12. 


2.41 x 10~ 8 


8.3 x 10~ 9 


3.6 x 10~ 9 


8. x 10" 10 






13. 


1.02 x 10~ 8 


4.4 x 10~ 9 


1.7 X 10 -9 


3. x 10" 10 






2.5 


4.15 x 10~ 3 


5.0 x 10" 4 


5.5 x 10" 4 


2.3 X 10" 4 






3.5 


4.54 x 10~ 4 


3.2 x 10~ 5 


4.9 x 10~ 5 


2.5 x 10~ 5 






4.5 


5.07 x 10~ 5 


4.8 x 10~ 6 


5.2 x 10~ 6 


2.8 x 10~ 6 






5.5 


1.14 x 10~ 5 


1.8 x 10~ 6 


1.1 x 10~ 6 


6. x 10" 7 


UJ 




6.5 


3.33 x 10~ 6 


4.5 x 10~ 7 


3.8 x 10~ 7 


1.8 x 10~ 7 






7.5 


7.7 x 10~ 7 


2.0 x 10~ 7 


1.0 x 10~ 7 


4. x 10" 8 






9. 


1.94 x 10~ 7 


6.4 x 10~ 8 


2.9 x 10~ 8 


1.1 x 10" 8 






11. 


4.8 x 10~ 8 


1.6 x 10~ 8 


1.0 x 10~ 8 


3. x 10" 9 
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TABLE XI: The invariant differential cross section 


— ^ — H d , a of Kg, r/, and meson production measured 


111 Lilt lllUlUdLt^U 


decay channel. Notations are: 


V is the differential cross section, 


A, B, and C are 


the three types of errors 


described in the 


text. 












meson decay 


Pt 


y 


A 


B 


C 


channel 


GcV 1 c 




mb/(GcV/c) 2 






2 25 


A fifi v i n — 3 

L t. DO X _LU 


7.7 x 10~ 4 


7.4 x 10~ 4 


3.0 x 10 -4 




2 75 


i v in -3 

l.OU A ±u 


1.1 x 10~ 4 


1.8 x 10~ 4 


8.0 x 10~ 5 




3 25 


4 14 X 10 — 4 


2.5 x 10~ 5 


5.5 x 10~ 5 


2.6 x 10 -5 




3 75 


1 54 x 10 — 4 


8.0 x 10~ 6 


2.0 x 10~ 5 


1.0 x 10~ 5 




4 25 


^ riQ v- i n — ^ 
o.uy a j.u 


2.8 x 10~ 6 


6.6 x 10~ 6 


3.2 x 10~ 6 




4 75 


9 99 v- in - 5 


1.2 x 10~ 6 


2.9 x 10~ 6 


1.4 x 10" 6 






l.UD X 1U 


6.0 x 10~ 7 


1.4 x 10~ 6 


7.0 x 10~ 7 




5.75 


4.74 X 10 — 6 


3.3 x 10~ 7 


6.1 x 10~ 7 


3.0 x 10" 7 




6.25 


2.74 x 10 — 6 


2.2 x 10~ 7 


3.6 x 10~ 7 


1.7 x 10~ 7 




6.75 


1.30 X 10 — 6 


1.3 x 10~ 7 


1.7 x 10~ 7 


8.0 x 10~ 8 




7.25 


7.70 X 10 — 7 


1.0 x 10~ 7 


1.0 x 10~ 7 


5.0 x 10" 8 




7.75 


3 82 -x 1 -7 


6.0 x 10~ 8 


5.3 x 10~ 8 


2.4 x 10~ 8 




8.25 


9 fiS y 1 n -7 


4.4 x 10~ 8 


4.1 x 10~ 8 


1.8 x 10" 8 




8.75 


1 59 x 1 -7 


3.1 x 10~ 8 


2.3 x 10~ 8 


1.0 x 10~ 8 




9 50 


7 an v in -8 

l.OU A 1U 


1.2 x 10~ 8 


1.1 x 10~ 8 


4.9 x 10~ 9 




10 5 


s 4Q v in -8 

O . A _LU 


8.5 x 10~ 9 


5.6 x 10~ 9 


2.2 x 10 -9 




11.5 


1.25 x 10 — 8 


3.7 x 10~ 9 


2.2 x 10~ 9 


8.0 x 10- 10 




12.75 


6 50 -X 1 -9 


2.1 x 10 -9 


1.3 X 10" 9 


4.1 x 10 _1 ° 




■~( 9^; 
o.zo 


l.Do X IU 


3.6 x 10" 5 


1.6 X 10" 5 


7. x 10" 6 




3.75 


7.38 X 10 — 5 


8.9 x 10~ 6 


7.3 x 10~ 6 


3.3 x 10~ 6 




4.25 


2.55 x 10~ 5 


4.1 x 10~ 6 


2.5 x 10~ 6 


1.1 x 10^ 6 




4.75 


i v in -5 

l,Ol7 A ± U 


2.0 x 10~ 6 


1.3 x 10~ 6 


6. x 10~ 7 




5 25 


4 9 x 10 — 6 


1.0 x 10~ 6 


5. x 10~ 7 


2. x 10 -7 




5.75 


9 Q9 y in -6 
Z, . O.^ A _LU 


4.3 x 10~ 7 


2.3 x 10~ 7 


1.0 x 10~ 7 


Tjf Tj7r~ >r 7T _ 


6.25 


1.13 x 10 -6 


3.1 x 10~ 7 


1.1 x 10~ 7 


5. x 10~ 8 




6.75 


7.7 x 10 -7 


1.7 x 10~ 7 


8. x 10~ 8 


3. x 10 -8 




7.5 


2.33 x 10 — 7 


5.4 x 10~ 8 


2.7 x 10~ 8 


1.0 x 10~ 8 




8 5 


1 07 x 10 — 7 


3.3 x 10~ 8 


1.2 x 10~ 8 


5. x 10 -9 




9.5 


5.2 x 10 -8 


1.7 x 10~ 8 


6. x 10~ 9 


2. x 10~ 9 




10.75 


1.09 x 10 -8 


5.8 x 10 -9 


2.1 x 10 -9 


4.9 x 10 _1 ° 




0.125 


0.264 


6.3 x 10 -2 


2.6 x 10 -2 


1.1 x 10 -2 




0.375 


0.188 


3.1 x 10~ 2 


1.8 x 10~ 2 


8. x 10~ 3 




0.625 


8.9 x 10 — 2 


1.5 x 10~ 2 


9. x 10~ 3 


4. x 10 -3 




0.875 


5 83 y 1 -2 

U.UJ A 1U 


8.2 x 10" 3 


5.8 x 10~ 3 


2.5 x 10~ 3 




1.125 


9 W y If) -2 


4.3 x 10 -3 


2.5 x 10 -3 


1.2 x 10 -3 


' + — 
Tj Tj7T 7T 


1.375 


1.31 x 10 — 2 


2.7 x 10 -3 


1.3 x 10 -3 


6. x 10" 4 




1.75 


2 79 x 10 — 3 


7.5 x 10~ 4 


2.8 x 10 -4 


1.3 x 10 -4 




2.5 


7.2 X 1U 


1.5 x 10~ 4 


7. x 10 -5 


3. x 10 -5 




3.5 


9.7 X 10" 5 


3.1 x 10" 5 


1.0 x 10" 5 


4. x 10" 6 




1.1 


3.32 x 10 -2 


2.6 x 10 -3 


2.6 x 10 -3 


4. x 10 -4 




1.45 


1.01 x 10~ 2 


7. x 10~ 4 


5. x 10~ 4 


1. x 10" 4 




1.95 


3.16 X 10~ 3 


1.9 x 10~ 4 


1.7 x 10 -4 


4. x 10~ 5 




2.45 


9.28 x 10~ 4 


5.6 x 10~ 5 


5.3 x 10~ 5 


1.1 x 10 _s 




2.95 


2.99 x 10~ 4 


1.9 x 10~ 5 


1.8 x 10~ 5 


4. x 10~ 6 


<p K + K~ 


3.45 


1.02 x 10~ 4 


6. x 10~ 6 


6. x 10 -6 


1. x 10~ 6 




3.95 


3.49 x 10~ 5 


2.6 x 10~ 6 


2.3 x 10 -6 


4. x 10" 7 




4.45 


1.38 x 10~ 5 


1.8 x 10~ 6 


9. x 10" 7 


2. x 10~ 7 




5.5 


2.31 x 10" 6 


4.1 x 10~ 7 


1.6 x 10" 7 


3. x 10" 8 




7. 


3.21 x 10~ 7 


7.9 x 10~ 8 


2.4 x 10~ 8 


4. x 10~ 9 
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